Metastability in Spin-Polarized Fermi Gases 
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We study the role of particle transport and evaporation on the phase separation of an ultracold, spin-polarized 
atomic Fermi gas. We show that the previously observed deformation of the superfluid paired core is a result 
of evaporative depolarization of the superfluid due to a combination of enhanced evaporation at the center of 
the trap and the inhibition of spin transport at the normal- superfluid phase boundary. These factors contribute 
to a nonequilibrium jump in the chemical potentials at the phase boundary. Once formed, the deformed state is 
highly metastable, persisting for times of up to 2 s. 
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The BCS theory of superconductivity is remarkably suc- 
cessful in describing pairing of unpolarized spin-72 parti- 
cles. Pairing in spin-polarized systems is much more compli- 
cated, however, prompting speculation about exotic new pair- 
ing mechanisms that began shortly after the development of 
the BCS theory [1, 2], and continues until today [3, 4]. Spin- 
polarization, or more generally, imbalanced Fermi energies, 
arise in several physical situations including certain supercon- 
ductors that support coexisting magnetic and superconducting 
order, color superconductivity in quark matter, and in ultra- 
cold atomic gases created with imbalanced spin populations. 
In 2006, a group at Massachusetts Institute of Technology 
(MIT) [5, 6] and our group at Rice University [7, 8] discov- 
ered that strongly-interacting spin-imbalanced trapped atomic 
gases undergo a first-order phase separation between a fully- 
paired superfluid core and lower density polarized regions. 

There are significant qualitative and quantitative differences 
between the MIT and Rice experiments. The phase separa- 
tion in the case of the Rice experiment was characterized by 
strong deformation of the paired core, in violation of the local- 
density approximation (LDA). In the LDA, the local chemical 
potentials jUo-(r) = /i, T - V(r) depend only on the trap potential 
V(r) and the spatially- uniform global chemical potentials fi tT . 
Here, <x =T,4 designate the two states of a pseudo-spin- 1 /2 
system. The local densities are given by these local chemi- 
cal potentials and the equilibrium equation of state of an in- 
finite, spatially uniform system. For a harmonically confined 
gas with an unpolarized central core, the LDA implies a flat- 
topped axial spin-density (obtained by integrating the three- 
dimensional spin density along the two radial coordinates) [9- 
11]. While this flattopped distribution was observed in the 
MIT experiment [6], in the Rice experiment the shape of the 
paired core was significantly less elongated than V(r), result- 
ing in a central dip in the axial spin-density [7, 8]. Further- 
more, pairing in the Rice experiment was much more robust 
than in the MIT experiment, persisting to much larger popu- 
lation imbalances. This robust pairing is apparently in con- 
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tradiction with the Clogston-Chandrasekhar limit describing 
the break-down of pairing when the difference between the 
chemical potentials of the two spin-states exceeds the pairing 
gap [12-15]. 

Possible explanations for this discrepancy have focused on 
the primary differences between the two experiments, which 
are trap aspect ratio and particle number [11, 16-21]. In both 
cases, the confining potential is approximately harmonic and 
elongated along the cylindrical (z) axis, as shown in Fig. 1(a). 
For the Rice experiment, however, the ratio of the radial to ax- 
ial trap frequencies was ~30, while for MIT it was ~5. Also, 
the total particle number in the Rice experiment was ~10 5 , 
while for MIT it was ~5 x 10 6 . It was shown that while the ob- 
served deformation was consistent with the effect of a strong 
surface tension at the superfluid-normal interface [8, 16, 17], 
the required magnitude of the surface tension was inconsis- 
tent with detailed calculations [22, 23]. Furthermore, a recent 
experiment is largely in agreement with the MIT results de- 
spite having an aspect ratio and particle number that are sim- 
ilar to the Rice experiment [24]. A new mechanism has been 
proposed [25] which has its origins in the inhibition of ther- 
mal [26] and spin transport [25] across the phase boundaries 
coupled with an enhanced probability for evaporation at the 
axial center of the trap. These factors enable a distribution 
that is out of chemical equilibrium, where the difference in 
chemical potentials is depressed in the superfluid phase at trap 
center relative to the polarized normal phases in the wings. 
We have performed several experiments that strongly support 
this conjecture. In addition, we find that, once produced, the 
deformed state is remarkably metastable. 

We produce imbalanced mixtures of 6 Li atoms as before in 
the two lowest energy hyperfine states (F = m F = 
and (F = l /2, m F = - l /2), designated as states |T/> and |J,), 
respectively [7, 8]. A bias magnetic field is tuned to 834 G, 
corresponding to the unitary limit of the 6 Li Feshbach reso- 
nance. The atoms are confined in a hybrid optical-magnetic 
trap formed by a single focused laser beam propagating along 
the direction of the bias field (axial direction). Radial con- 
finement is produced by the Gaussian intensity profile of the 
laser beam, while axial confinement arises from the combi- 
nation of the Lorentzian axial profile of the laser beam and 
the residual magnetic curvature (confining) from the slightly 
non-Helmholtz configuration of the magnetic bias coils. This 
combination, depicted in Fig. 1(a), results in an aspect ratio 
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FIG. 1. (Color online) Plots of the hybrid magnetic-optical trap. 
Atoms are trapped at the focus of a laser beam propagating in the z 
(axial) direction. The laser wavelength is 1 .08 /jm and the beam is fo- 
cused with a Gaussian beam waist ( 1 /e 2 radius) of 30 fim. A residual 
magnetic curvature contributes to the axial confinement with a har- 
monic frequency of 3.8 Hz. (a)-(c) Optical trap depth U = 1.68/iK, 
corresponding to t = 700 ms for the 1 s evaporation trajectory. U g is 
the effective trap depth accounting for gravity. The relative contribu- 
tion of the magnetic curvature and gravity are small, (d)-(f) Potential 
at the final optical trap depth Uf = 0.74 fiK at t = 1 s, where the rel- 
atively strong magnetic curvature has the effect of opening up a lip 
at z = 0. At this trap depth, the combined axial frequency, due to the 
optical and magnetic forces is 4.7 Hz. Values of k B T F from Fig. 3(a) 
are indicated by the horizontal blue lines. 



of 90~100 with an approximately isotropic trap depth when 
the optical potential dominates. Evaporative cooling is ef- 
fected by lowering the laser power, such that the trap depth, 
as well as the trap aspect ratio, is gradually reduced. At suf- 
ficiently low optical power the magnetic curvature dominates 
the axial confinement. When this happens, the trap depth be- 
comes anisotropic, with the depth being largest along the axial 
(magnetic curvature) direction. Thus, a "lip" of minimum trap 
depth is formed in the radial direction at the axial trap center 
(z = 0), as shown in Fig. 1(d). Furthermore, gravity reduces 
the trap depth in the direction pointing downward in the lab 
(along -y). 

Sequential absorption images [8] record the column den- 
sity distributions of the trapped atoms for each state. Fig- 
ure 2 shows representative axial density profiles, obtained 
by integrating the column density images along the remain- 
ing radial coordinate, for images recorded at various times 
along the evaporation trajectory used in our previous studies. 
For this trajectory, the trap depth was reduced exponentially 
as e~' lr from its initial value I/, = to a final value 

Uf = 0.74yi/K in a total time t lot — 1 s, and with an expo- 
nential time constant t = 200 ms. Deformation is evidenced 
by a dip in the axial spin density, which begins to develop at 
approximately t = 700 ms [Fig. 2(b)], corresponding to a tem- 
perature T ^ 0.2 Tf, where Tf is the Fermi temperature of 
a noninteracting trapped gas of |f) atoms. The deformation, 
characterized by the parameter a = n v jn p [see Fig. 2(c)], in- 



creases as the evaporation progresses and is maximum at t tot , 
where T ^ 0.06 7> is at its minimum. 

Figure 3 shows the progression of several relevant param- 
eters during evaporation. U g , shown in Fig. 3(a), is the trap 
depth including gravity, as defined in Fig. 1(c). Also shown 
in Fig. 3 are parameters extracted from a less aggressive evap- 
oration trajectory for which [/, = 160 fiK as before, but now 
with Uf - 2.2 fiK, t to t - 3.4 s, and r = 500 ms. This trajectory 
is designed to be similar to the final part of the trajectory used 
in Ref. [24], where no deformation was observed. Figure 3(b) 
shows the value of rj — {U g - e p )/ksT, where e p = l /2tna£R 2 , 
u) z is the axial trap frequency, and R is the axial radius where 
the density of the majority state (f) goes to zero. The value of 
T] is an approximate measure of the closeness of the chemical 
potential to the trap lip, and hence is related to the rate of evap- 
oration. This quantity is significantly larger for the "gentle" 
3.4 s trajectory as compared to the "aggressive" 1 s trajectory, 
indicating a much slower rate of evaporation. Nonetheless, 
even though the 3.4 s trajectory is not as deep or as aggressive, 
the final temperature of ~0.09 Tp is similar to that achieved 
with the 1 s trajectory. Furthermore, as shown in Fig. 4, the 
axial spin density at the end of the trajectory is flattopped, in- 
dicating that there is no deformation even though the trap as- 
pect ratio at the end of evaporation is highly elongated (aspect 
ratio of ~96). Deformation is prevented in the 3.4 s trajectory 
by its higher final trap depth, which both reduces the rate of 
evaporation and minimizes the lip at z = 0. 

To determine whether the deformed state is only dynami- 
cally stable, existing only during rapid anisotropic evapora- 
tion, or rather is a metastable state, we ramped the trap depth 
up over a time period of 600 ms following evaporation, as 
shown in Fig. 5(a). This serves to significantly suppress the 
rate of evaporation, as can be seen from the nearly constant 
value of Tf in Fig. 5(a) and the large values of r\ in Fig. 5(b). 
Nonetheless, Fig. 5(d) shows that the deformation a remains 
for more than 2 s following trap recompression. The degree 
of deformation is seen to decrease following recompression, 
roughly on the same timescale of an observed rise in the tem- 




1000 500 1000 500 1000 500 

z(nm) 

FIG. 2. (Color online) Axial densities taken at various times t along 
the evaporation trajectory. The global polarizations P are: (a) 600 
ms, P = 0.14, (b) 700 ms, P = 0.21, (c) 740 ms, P = 0.28, (d) 1 
s, P = 0.18. The variation in P is a result of shot-to-shot variations, 
as each image requires the trap to be reloaded and evaporated to the 
specified /. The upper (black) curves correspond to the majority state 
(It)), the middle (blue) curves to the minority state (||», and the axial 
spin densities (|f > - ||» are given by the lower (red) curves. 
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FIG. 3. (Color online) Parameters extracted from the axial densi- 
ties for both the (aggressive) 1 s and (gentle) 3.4 s trajectories. T is 
determined from the mean energy E of unpolarized (P = 0) distri- 
butions which are separately evaporated using the same trajectory. E 
for each distribution is obtained from its mean-squared radius via the 
virial theorem [27]. The E vs T calibration is given in Ref. [28] and 
is based on the experimental data of Ref. [29]. 77 = (U g - e p )/k B T, 
where e p = l l2mu>zR 2 , a> z is the axial trap frequency, and R is the 
axial radius where the density of the majority state (T) goes to zero. 
The deformation parameter is defined as a = n v /n p , as depicted in 
Fig. 2(c). The error bars are mainly statistical uncertainty from the 
average of ~6 shots of various values of P at each value of /. The 
vertical dashed line in (d) indicates the onset of deformation. 



perature. 

Another measure of the effect of the lip in the potential may 
be obtained by axially displacing the center of the magnetic 
curvature with respect to the focus of the optical trap laser 
beam, as depicted in Fig. 6. Since the lip is located at the 
minimum of the magnetic curvature (z = z,„), its position no 
longer coincides with the overall minimum of the combined 
magnetic-optical potential, indicated by zo in Fig. 6. Figure 6 
shows that while the unpaired atoms, given by the distribution 
°f IT) _ reside near the center of the combined potential 
(zo), the paired core, given by the \i) distribution, displaces 
towards the lip at z = z m where evaporative cooling has max- 
imum effect. This observation provides a graphic illustration 
of the lack of equilibration between the superfiuid core and 




FIG. 4. (Color online) Axial densities for the 3.4 s trajectory at t = 
3.4 s, with P = 0.24. Curve designations are the same as in Fig. 2. 
The dashed vertical lines indicate the location of phase boundaries. 
The flattopped axial spin density is consistent with the LDA, even 
though the aspect ratio of the trap potential is ~96. 
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FIG. 5. Recompression of the trap following the 1 s evaporation. 
Here, the 1 s evaporation finishes at t = and is followed by a slow 
recompression of the trap over the next 0.6 s. Even though the fi- 
nal value of U g is similar to that shown for the 3.4 s trajectory in 
Fig. 3 the recompressed values of T F are lower due to smaller overall 
numbers (~3 x 10 4 vs ~10 5 ). 
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FIG. 6. (Color online) Column density images for an axial displace- 
ment of the center of magnetic curvature located at z m = and in- 
dicated by the dotted vertical line, from the focus of the optical trap 
beam located at z = -1 mm. At this trap depth (U = 0.65 /uK), 
the center of the combined potential, indicated by the dashed vertical 
line, is located at zo = -210^m. The trap beam waist is 26 yum and 
the residual magnetic curvature is 6.4 Hz for these data. The image 
size is 1653 fim x 100 /jm and P = 0.63. The uncertainty in zo and z„, 
is 20 /im due to uncertainties in the optical trap parameters. 



the normal phases. 

We have presented a series of measurements that are con- 
sistent with a model of evaporative depolarization [25]. In 
this model, the chemical potential of the majority species is 
depleted by preferential evaporation in the vicinity of a lip in 
the elongated trap potential. Because of the inhibition of spin 
transport, for sufficiently fast evaporation the chemical poten- 



tial difference can be smaller in the paired core than in the 
polarized wings. This helps to stabilize the superfiuid phase 
beyond the Clogston-Chandrasekhar limit, as was observed 
experimentally [7, 8]. The absence of a spatial variation in the 
spin-density of the superfiuid phase [8] shows that the gradient 
in the chemical potential difference is not smoothly varying, 
but rather that it undergoes a sudden jump at the phase bound- 
ary. This indicates that the interface, not the bulk superfiuid, 
is the dominant obstacle to spin transport. Slow relaxation of 
a nonequilibrium spin distribution in a spin-imbalanced phase 
separated gas was also reported in Ref . [30] . The observed re- 
laxation time scales in that experiment are consistent with ours 
but the relative importance of the interface compared with the 
bulk superfiuid could not be determined. (A previous exper- 
iment also reported slow diffusion but it was conducted with 
a balanced spin mixture, and hence without phase boundaries 
[31]). In our experiment, the jump in the local chemical po- 
tential difference at the phase boundary strongly deforms the 
shape of the superfluid-normal interface, making it much less 
curved than it would be at equilibrium. This nonequilibrium 
distribution is remarkably metastable, consistent with the cal- 
culations presented in Ref. [21], in which they find both LDA- 
like solutions as well as nearly degenerate LDA-violating ones 
whose density distributions closely resemble those that we ob- 
serve. 
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